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Abstract 

The CDF experiment reported a lepton photon missing transverse energy {^t) 
signal 3a in excess of the Standard Model prediction in Tevatron Run I data. The 
excess can be explained by the resonant production of a smuon, which subsequently 
decays to a muon, a photon and a gravitino. Here, we perform combined fits of this 
model to the CDF ^l^x excess, the DO measurement of the same channel and the 
CDF 7-^7^ channel. Although the rates of the latter two analyses are in agreement 
with the Standard Model prediction, our model is in good agreement with these data 
because their signal to background efficiency is low at the best-fit point. However, 
they help to constrain the model away from the best fit point. 



1 Introduction 

The CDF experiment has recently discovered an anomaly in the production 
rate of lepton-photon-^^-r in pp collisions at y/s = 1.8 TeV, using 86.34 pb~^ 
of Tevatron 1994-95 data [1]. While the number of events expected from the 
Standard Model (SM) were 7.6 ± 0.7, the experimentally measured number 
corresponded to 16. Moreover, 11 of these events involved muons (with 4.2 ± 
0.5 expected) and 5 electrons (with 3.4 ± 0.3 expected). 

In earlier papers [2] we suggested that the excess can be simply understood in 
terms of the minimal supersymmetric standard model (MSSM) which has the 
following ingredients: (1) the model is i?-parity violating with an L-violating 
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Fig. 1. Feynman diagram of resonant smuon production followed by neutralino de- 
cay. 

X211 coupling, and (2) the supersymmetric spectrum includes an ultra-light 
gravitino of mass ~ 10~^ eV. We have demonstrated in our earlier papers [2] 
that such a model provides a natural explanation of the CDF anomaly and 
explains not only the excess in the number of H'J^t events but also explains 
the main features of the kinematic distributions of the excess events. The 
excess can be explained using a small value of the L-violating coupling A211 ~ 
0.01 because of the resonant production of smuons in the annihilation of an 
initial-state qq pair. The smuon thus produced decays predominantly into a 
bino-dominated neutralino and a muon, with the neutralino further decaying 
into a photon and a gravitino. The production and decay has been shown 
in the Feynman diagram in Fig. 1. The fact that the excess is seen in final 
states involving photons emerges very neatly in our model because the decay 
Xi iG dominates overwhelmingly over other decay modes. 

In the current article, we extend the previous studies by including two addi- 
tional pieces of independent empirical information. We include the DO Run I 
measurement [3] of the /i7 missing Et process, as well as 7 missing Et data 
coming from CDF in Run I [4]. The empirical background event-rate in the 
fj-j^T channel is quite different in the CDF and DO cases due to the different 
cuts employed. Our scenario predicts excesses in each of these channels, and 
we determine to what extent it is in accord with their measured rates. By 
performing a combined fit to all three event rates, we constrain the masses of 
the relevant sparticles in the event, as well as \'2ii- 



2 The model 

In order that the cross-section for the production of the smuon resonance is 
substantial enough to account for the anomalous events, we need the left- 
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Fig. 2. Feynman diagram of resonant muon sneutrino production followed by neu- 
tralino decay, resulting in a signature. 

handed smuon to be light i.e. m^^ ~ 150 GeV. Further, one needs to couple 
the smuon to valence quarks in the initial state implying that the L- violating 
operator that we need is L2Q1D1 corresponding to the coupling A211, which 
generates the interactions jlud and z/^dd (and charge conjugates), along with 
other supersymmetrised copies involving squarks. Therefore, the operator we 
invoke in our model predicts supersymmetric signals in other channels which 
manifest themselves through the production of either sneutrinos or squarks. In 
our model, we take the squarks to be heavy and so their effects on experimental 
observables will be negligible. On the other hand, the spontaneously broken 
SU(2)l symmetry in the MSSM implies that muon sneutrinos have a tree-level 
mass squared [5] 

'T^l = - "^J + (1 - sin^ e^) cos(2/3)M| tjh, < m^^ , (1) 

where is the mass of the muon, the Weinberg angle, and tan j3 the ratio 
of the MSSM Higgs vacuum expectation values (cos 2/3 < 0). Since < m^^ 
if we assumed that smuons were produced at the Tevatron Run I energy, we 
should expect that muon sneutrinos could also be produced. The dominant 
production mechanism is resonant sneutrino production and subsequent decay, 
as shown in Fig. 2. It results in a experimental signature. 

The pattern of masses of the super-particles suggested by the CDF fi'y^j' 
anomaly is the following: the smuon is around 150 GeV, and the only other 
light sparticles are the neutralino (which is typically about 45 GeV lighter than 
the smuon) and the ultra-light gravitino (which is as light as 10"'^ eV). We 
enforce degeneracy between the first two generations in order to avoid flavour 
changing neutral currents. Other sparticles do not play a role in this analysis, 
and are set to be arbitrarily heavy. Such a light gravitino materialises naturally 
in models of gauge-mediated supersymmetry breaking (GMSB)[6]. However, 
in the minimal version of GMSB models the chargino and the second-lightest 
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neutralino are not much heavier than the neutrahno and this feature of the 
minimal model is certainly not desirable for our considerations because it leads 
to large jets+7+^y rates which are not seen by experiments ^ . If the minimal 
versions of GMSB models do not yield the pattern of super-particle masses 
that we need, then the obvious question to ask is what is the high-energy 
model that yields this mass spectrum at low energies. It is interesting to note 
that such a mass spectrum can arise in an alternate model of GMSB which 
is obtained from compactifying 1 1-dimensional M-theory on a 7- manifold of 
G2 holonomy [10]. For the purposes of this paper, however, we simply work 
with the low-energy model with the mass spectrum described above and do 
not worry about the high energy completion of this model. 



3 Simulating the experiments 

In our model, we have essentially four free parameters that are relevant to the 
data we fit: the gravitino mass, ttiq, the neutralino mass m^o, the smuon mass, 
Mji and the i?- violating coupling Xi^n- However, instead of simultaneously 
fitting the four parameters using the experimental data, we choose to work 
with fixed values of two of these parameters close to their best-fit values while 
performing fits in the two other parameters. For our analysis, we take other 
parameters like tan/? = 10 to be constant. The coupling, Agn, is constrained 
fi:om = r(7r eu)/{'K ixv) [11] to be < 0.059 x [12]. But since 

the constraint involves a squark mass which is large in our model, it is not 
very relevant. While the production of the smuon resonance is through the 
/^-violating mode, to produce the /7^t excess we require that its decay goes 
through the i?-conserving channel to a neutralino and muon final state. The 
it!-violating decay of the slepton is possible but constrained, in principle, by 
the Tevatron di-jet data [13] which exclude a a.B > 1.3 x 10^ pb at 95% C.L. 
for a resonance mass of 200 GeV. However, in practice this does not provide 
a restrictive bound upon our scenario as long as the i?-violating coupling is 
sufficiently small, < 0{1). We also add that the di-jet bound is not restrictive 
because it suffers from a huge QCD background. By restricting A211 to be 
small, we also avoid significant rates for the possible i?- violating decays of 
Xi or X? ^n- 

We use the ISASUSY part of the ISAJET7 . 58 package [14] to generate the spec- 
trum, branching ratios and decays of the sparticles. For an example of param- 
eters, we choose (in the notation used by ref. [14]) tan/5 = 10 and A-t^T^f, — O.fi 
together with other fiavour diagonal soft supersymmetry breaking parameters 



^ This observation has been made earlier in the Uterature [7] in the context of the 
GMSB-based explanation [8] of the eej-f^T CDF event [9]. 
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£:T(7)/GeV 


<60 


60-65 


65-80 


>80 


efficiency 


40% 


47% 


51% 


54% 



Table 1 

Pfioton efficiency in the CDF 'y^j- analysis 



experiment 


luminosity 


observed number 


background 


CDF 


86 pb-i 


11 


4.2 ±0.5 


DO iij^T 


73 pb-i 


58 


58 ± 9.75 


CDF 


86 pb-i 


11 


11 ±2.2 



Table 2 

Observed number of events passing cuts in the text and Standard Model back- 
grounds for the three pieces of data included in the combined fit. The integrated 
luminosity for each analysis is also listed. 

are set to be very large. We emphasise that this is a representative point in 
the supersymmetric parameter space and not a special choice. 

As stated earlier, we present analyses for three sets of data in this paper: 

the CDF Run I data on Ij^j- [1], the DO Run I data on fij^^ [3] and the 
CDF Run I data on 7-^2- [4]. We now present our fiducial efficiencies and cuts, 
which mimic those of the relevant experiments. The CDF experiment detects 
photons with the constraints that the following pseudo-rapidity regions are 
excluded: |77^| > 1, |r/^| < 0.05 and the region 0.77 < rj'^ < 1.0, 75° < (f)< 90°. 
For the Ij^t data, photon detection efficiency within these cuts is r)^ — 81%. 
Muons have a 60% detection efficiency if \ri^\ < 0.6 or 45% if 0.6 < f]^ < 1.1. 
To improve signal to background ratio efficiency, the following cuts are used: 
> 25 GeV, ^^(7) > 25 GeV and > 25 GeV. For the 7^^ analysis, 
the CDF experiment [4] has chosen the following cuts: \7]'^\ < 1.0, E]^ > 55 
GeV and > 45 GeV. The cut on the photon Et is chosen to be as large 
as 55 GeV so as to beat down the background due to cosmic rays. Photon 
fiducial efficiencies are shown in table 1, and were obtained from ref. [4] 

In contrast, in order to simulate the DO experiment, we specify < 1.0, 
\rij\ < 1.1 or 1.5 < jr^-yl < 2.5. We use the same cuts as DO: pf^ > 15 GeV, 
> 10 GeV on the muon and photon Pt respectively. Also, > 15 GeV, 
AR{fi-f) > 0.7 and Mrifi^r) > 30 GeV. Here, Ai?(/x7) is the distance be- 
tween the photon and muon in pseudo-rapidity {1]) and transverse angle {(f)) 
space. M|. = E^—p'^ is the transverse mass. Within these cuts, we have fidu- 
cial efficiencies of 71.1%, 50.1% and 51.0% for the trigger, muon and photon 
respectively. These cuts have been optimised by DO with a view to studying 
the effects of anomalous gauge boson couplings on this final state. Unfortu- 
nately, the signal to background ratio for our model in the DO analysis is far 
from optimal for the signal that we propose to study. 
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The observed number of events in the different analyses, and their Standard 
Model backgrounds (taken from the experimental papers [1,3,4]) are shown in 
Table 2. 

We now simulate the signal events for the processes in Figs. 1,2. We use 
HERWIG6.4 [15] including parton showering (but not including jet isolation 
cuts) to calculate cross-sections for single slepton production. 



4 Combined fits 

For each of the three data listed in Table 2, we can define a log likelihood 
defined by the Poissonian log-likelihood for Uo observed events, Us expected 
signal events (for fixed values of SUSY breaking parameters) , convoluted with 
a Gaussian probability distribution of the number of expected Standard Model 
background events usm and its uncertainty asu' 

ln£(ns) = — / e '^"sm (noln(n -|- ng) — {n-\- ng) — Inrio) dn.{2) 

This is a good approximation provided nsu is several times asu above zero, 
as is the case here. In(i2(ns)/i2(0)) is then the (signalH-background) to back- 
ground hkelihood ratio for a given analysis. We can form the total fit hkelihood 
ratio In by adding the likelihood ratio from each of the three analyses in 
Table 2. We will always consider two relevant parameters to be fixed and 
fit the model to the other two. The best-fit point in parameter space corre- 
sponds to the maximum value of the likelihood ratio and is found by using 
MINUIT [16]. — 21n>Cj^^3. corresponds to the equivalent number of (Standard 
Model) - (best-fit )= Ax^. We have one degree of freedom, and therefore 
the 90% and 95% confidence level (CL.) limits on parameters then lead 
to [16] ln£*°^^ — In/^*"* = 1.35, 1.92 respectively in the fit. We will constrain 
M^o > 100 GeV, as implied by jj^t LEP2 data [17]. We now discuss the 
results of combined fits for various hyper planes of parameter space. 

Fig. 3 displays the 90 and 95% C.L. fit regions as the area between the solid 
lines in the \2ii, M^o plane. Am = Mji — M^o has been kept fixed at 45 GeV in 
order to keep the decay mode p, —>■ nxi open and ttIq — 10~^ eV. A significant 
amount of parameter space fits the combined data, with ranges A211 > 0.001. 
Increasing M^o produces a lower cross-section because of kinematical suppres- 
sion, but this effect can be compensated by raising A211, thus increasing the 
production rate. The best-fit point is X211 = 0.11, M^o = 100 GeV, with 
Ax^ = 6.90. Fig. 3a shows that we expect between 3-8 signal events in the 
CDF I'f^x channel. This data dominates the fit because backgrounds (and 
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M^o (GeV) 

Fig. 3. Predicted number of excess events in (a) CDF Ij^j-, (b) DO l^^xi ^^^d (c) 
CDF 7-^^ channels in the A2ii-neutrahno mass plane for tan/3 = 10, Am = 45 GeV 
and uIq = 10^'^. Labelled contours of equal numbers of signal events are shown as 
dashed curves. The 90% C.L. and 95% C.L. combined-fit regions are between the 
inner and outer pair of solid lines respectively. 

their uncertainties) are larger in the other analyses. 1.5-6 events are expected 
at DO, and between 1.5 and 8 signal events were predicted for the CDF j^j, 
signature, depending on the parameters. 

Fig. 4 displays the 90 and 95% C.L. fit regions as the area to the left-hand side 
of the solid lines in the log(mg),M^o plane. Here, we use the default values 
A211 = 0.01, Am = 45 GeV. It is clear from the fit regions that the data prefer 
lower values of the neutralino mass. The best-fit point is: M^i = 100 GeV, 
rriQ = lO^'^ '^ eV, Ax^ = 6.76. The figure illustrates that if rriQ > 0.1 eV, 
the branching ratio of Xi ^ G'~f becomes tiny, decreasing the cross-section for 
the CDF I'y^T signal, which dominates the fit. When tuq is below 10~^ eV, 
fl — > fiG decays dominate, again decreasing the CDF fij^j- signal. 2-6 CDF 
Ij^T excess events, 2-4 DO excess events and 0-3 j^rp CDF excess events 
are expected, shown by the dashed contours within the 90% C.L. regions of 
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Fig. 4. Predicted number of excess events in (a) CDF l^^j-, (b) DO l^^f-, 
and (c) CDF 'y^rp channels in the gravitino mass-neutrahno mass plane, for 
tan (3 = 10, Am = 45 GeV and A211 = 0.01. Labelled contours of equal numbers of 
signal events are shown as dashed curves. The 90% C.L. and 95% C.L. combined-fit 
regions are to the left of the inner and outer solid lines respectively. 

Figs. 4(a-c). 

Fig. 5 displays the 90 and 95% C.L. fit regions as the area between the solid 
lines in the \'2ii, Mp^ plane. Here we have chosen default values of tUq = 10~^ 
100 GeV. A significant amount of parameter space fits the 



eV and M^o 



combined data. The best-fit point is 



211 



0.114, Mr, 



154 GeV, with 



Ax^ = 6.90. Fig. 5a shows that we expect between 1-12 signal events in the 
CDF I'j^j- channel. 1-8 events are expected at DO, and up to 4 signal events 
were predicted for the CDF ^^j, signature, depending on the parameters. 



Fig. 6 displays the 90 and 95% C.L. fit regions as the area enclosed by the solid 
lines in the m^, plane. Here we have chosen default values of X211 = 0.01 
and M^o = 100 GeV. The ranges = 130 - 230 GeV, = lO'^'^ - lO^^-^ 
eV provide a reasonable combined fit. The best-fit point is = lO^'^'^ eV, 




' 
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Fig. 5. Predicted number of excess events in (a) CDF Ij^j^, (b) DO I^^t^ ^^'^ (c) 
CDF 7-^2" channels in the A2ii-smuon mass plane, for tan (5 = 10, ttIq = 10~^ eV 
and M^o = 100 GeV. Labelled contours of equal numbers of signal events are shown 
as dashed curves. The 90% C.L. and 95% C.L. combined-fit region is between the 
inner and outer pair of solid lines respectively. 

Mfi = 142 GeV, with Ax^ = 6.76. Fig. 6a shows that we expect between 1-6 
signal events in the CDF Ij^q. channel. 0-4 events are expected at DO, and 
up to 0.2 signal events were predicted for the CDF signature, depending 
on the parameters. 



5 Conclusions 

We have provided combined fits for a supersymmetric model that explains the 
I'y^T CDF Run I excess in events, which was at the 2.7a level [1]. We have used 
the Run I 7.^^ data recently presented by CDF, as well as anomalous trilinear 
gauge boson coupling data from DO. Constraints upon various hyper planes in 
the gravitino, smuon, neutralino and R-parity violating coupling space have 
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Fig. 6. Predicted number of excess events in (a) CDF Ij^x^ 0^) -DO ^T-^Tj 
(c) CDF channels in the gravitino mass-smuon mass plane, for tan (3 = 10, 
M^o = 100 GeV and A211 = 0.01. Labelled contours of equal numbers of signal 
events are shown as dashed curves. The 90% C.L. and 95% C.L. combined-fit regions 
are depicted by solid lines. 

been displayed. In totality, the signal rates predicted by our model for the 
three analyses fit the data well, best fit points corresponding to a Ax^ = 6.9 
fit compared to the Standard Model, for one degree of freedom. 

Unfortunately, background rates in the DO anomalous trilinear gauge boson 
coupling data are too high for it to be very sensitive to the predicted signal 
rate. We note, however, that another analysis on existing DO Run I data with 
cuts optimised to test a ^■j^t excess would provide a good test of our scenario. 

The CDFj^rp channel suffers from a high -Et(7) > 55 cut because of cosmic 
backgrounds, which unfortunately also cuts most of the signal. It was shown 
in ref. [2], that if the cut could be reduced to 25 GeV, a signal rate of several 
times that in the CDF channel is possible. This is an important observa- 
tion for Run II, where cosmic backgrounds could be cut by additional timing 
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information in the detector. We look forward to the analysis of Run II data, 
which will be the final arbiter on l^^^ excess, as well as on our scenario. 
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